Observation and Characterization of the Vestige of the Jamming Transition in a 

Thermal 3D System 
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We study the dependence of the pair-correlation function, g(r), and particle mobility on packing 
fraction in a dense three-dimensional packing of soft colloids made of poly N-isopropyl acrylamide 
(pNIPAM), a thermo-sensitive gel. We find that g(r) for our samples is qualitatively like that of 
a liquid at all packing fractions. There is a peak in g±, the height of the first peak of g(r), as a 
function of packing fraction. This peak is identified as the thermal remnant of the T = divergence 
found at the jamming transition in simulations of soft frictionless spheres at zero-temperature. Near 
where there is a peak in g± the particles become arrested on the time scale of the experiment. 
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I. INTRODUCTION 

From hard-packed granular road beds to glass look- 
out ledges in skyscrapers, amorphous materials are used 
as solid support structures. Yet we do not understand 
the origin of rigidity in those materials. In the case of 
crystalline solids we understand the onset of rigidity as 
being a consequence of breaking translational symmetries 
during the fluid-to-crystal transition In contrast, no 
obvious symmetries are broken in the transition of a fluid 
to an amorphous solid. Is there any structural signature, 
however subtle, that can be identified with the onset of 
rigidity? Previous work has searched for such a signature 
in the pair-correlation function, g(r). Different criteria 
have been proposed. These include the splitting of the 
second peak of g(r) into two sub-peaks |2j,the ratio of 
the first minimum to the first maximum of g(r) [|], and 
changes to the contact-force distribution which is related 
to g(r) at small r |4(. 

Jamming is a way of constructing rigid amorphous 
structures [5j. Recent simulations of jamming at zero 
temperature [(j have suggested a different structural sig- 
nature for rigidity onset. The rigidification of random 
packings of finite-range soft repulsive spheres at tem- 
perature T = 0, is controlled by the packing fraction, 
(j). Such packings have a sharp jamming transition at a 
critical packing fraction, cf> c , where the particles first un- 
avoidably make contact. Below </> c no particles overlap, 
whereas above (j) c particles must overlap with their near- 
est neighbors and the system can support stress. This 
onset of rigidity is a purely geometric effect. At (j) c all 
nearest-neighbor pairs are separated by precisely one par- 
ticle diameter. This leads to a S- function in g(r) at its 
first peak. This divergence is a signature of the athermal 
(i.e., T = 0) jamming transition and varying above or 
below (j) c suppresses the divergence [6|. 

The relevance of jamming to systems at T > as well 
as its relation to dynamical arrest and the glass transition 
have not been clear. Dense suspensions of hard colloids 
have been thoroughly studied in the context of the col- 



loidal glass transition and super-cooled fluids 0, 
However, hard colloids are ill suited to study many as- 
pects of jamming since they cannot overlap so can only 
access configurations on the fluid side of jamming, that is 
(j) < cj) c . In this paper, we study as a function of <p a three- 
dimensional packing of soft colloids made from a thermo- 
sensitive hydrogel. Our particles are small enough to un- 
dergo Brownian motion and thus can be considered to be 
at finite temperature. To search for a structural signa- 
ture of the jamming transition, we use optical confocal 
microscopy and particle tracking techniques to determine 
the positions and displacements of particles in order to 
calculate the pair-correlation function, g(r), and the par- 
ticle mobility. 

In order to investigate the structure of our samples 
as a function of packing fraction we control <p by vary- 
ing the particle diameter at a fixed number density. 
As we vary 0, the height of the first peak of g(r), gi, 
varies in a non-monotonic fashion and has a maximum 
at (j) = 0*, near where there is also a dramatic reduction 
in the particle mobility. This maximum is a vestige of 
the zero-temperature jamming transition. The value of 
</>* is greater than C , the packing fraction where jam- 
ming would occur at T = 0. These observations are 
qualitatively consistent with recent simulations [19| and 
two-dimensional experiments using bi-disperse soft col- 
loids m. 

The fact that the first peak in g(r) is largest at </>* 
suggests that the sample at that <p may have more long 
range order than at other packing fractions. This can be 
checked by looking at the subsequent peaks of g(r) which 
contain information about medium- and long- range or- 
der. Close to and above 0* we observe that g(r) has at 
least 14 evenly spaced peaks which decay in height as r 
increases. In this regime the subsequent peaks do not 
vary appreciably with <\>. However, lowering <p below </>* 
strongly increases the damping of the higher-order peaks. 
We see no evidence of a split second peak in g(r) as is 
seen in hard-colloidal systems at the colloidal glass tran- 
sition (HI. 
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FIG. 1. a) The hydrodynamic diameter of the pNIPAM col- 
loids used in the experiments. Dashed line is an empirical fit 
to the data (solid circles). The horizontal line is the average 
nearest neighbor spacing for a given number density, n. b) 
Relative packing fraction. Vertical line in both graphs is the 
temperature where </> = </>* 



We also compare the structure and mobility of our 
samples as a function of packing fraction. We find a 
pair of systems on opposite sides of 0* with an order 
of magnitude difference in the particle mobility where 
the pair-correlation functions are experimentally indis- 
tinguishable. This suggests that the dynamics cannot be 
predicted simply from a knowledge of the structure. 



II. METHODS /MATERIALS 

We synthesize poly N-isopropyl acrylamide (pNI- 
PAM) colloids using surfactant-free emulsion polymer- 
ization [22], [23|| . The size, stiffness, and temperature de- 
pendence of the colloids is highly sensitive to the details 
of the synthesis and each batch must be calibrated in- 
dependently. The hydrodynamic particle diameter, cr, 
depends on the environmental temperature, T. We mea- 
sure <j{T) by observing the diffusion of the colloids in a 
dilute sample. As shown in Figure QJi, over the range 
24°C< T < 34°C, the diameter varies by nearly a factor 
of two. The data for cr(T) is empirically fit by two linear 
segments connected by a cubic spline. Above 35°C, the 
colloids collapse to a constant size[22|. The sample tem- 
perature, measured to a precision of 0.1°C, is controlled 
with a Bioptechs objective heater thermally coupled to 
the sample via the objective immersion fluid. 

In this paper, we estimate the packing fraction, 0, of 
dense systems from the particle diameter, cr and the num- 



FIG. 2. The corrected and raw g(r) curves for a sample at 



ber density, n, 



<XT) = -na(Tf. 



(1) 



Over a 10°C temperature range, <fi varies by a factor of 7 
as shown in Figur^Ib- We take data in the range between 
T = 27°C to 32°C. Over this range,the sample transforms 
from a liquid state, where the particles are diffusive, to 
a fully arrested packing. The primary effect of changing 
the environmental temperature is to change the size of 
the colloids (and hence <p) rather than the thermal energy. 

Sample chambers are made from microscope slides and 
cover slips and sealed with epoxy (Norland 61). The 
chambers are small, 0.15mm x 5mm x 3mm, to minimize 
thermal gradients within the sample. To prepare dense 
samples, sparse suspensions of colloids are centrifuged 
to sediment the colloids. The suspension is then heated 
to 35°C to de-swell the particles. The sediment is then 
pipetted into the sample chamber which is quickly sealed 
with epoxy so that the total volume, particle number, 
and number density, n, are fixed. After enclosure in the 
chamber, the samples are prepared for observation by 
first heating to 35°C to turn the system into a fluid. The 
sample is then cooled to 4°C in less than 5s by placing 
the sample chamber on 4°C metal surface. This quench 
protocol consistently generates arrested amorphous con- 
figurations. 

Particle dynamics and g(r) are extracted from data 
sets taken at different values of <p. Between subsequent 
measurements, the sample is melted and re-quenched. 
The waiting time from the quench to observation is ap- 
proximately 900 seconds and is controlled to minimize 
possible complications due to aging pll [25|. In order to 
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sample the (j) dependence of g(r) in finer increments the 
packing fraction was also swept continuously. For in- 
creasing (f) ramps, a fluidized sample is placed on a pre- 
heated objective. The objective heater is then turned off 
and the sample is periodically imaged as it cools slowly 
to ambient temperature. 

We image a two-dimensional slab far from the bound- 
aries of a three-dimensional sample using a Yokogawa 
CSU-XI confocal head and a Nikon 60x water immersion 
objective. Images are acquired at frame rates between 
0.3Hz and 10Hz using a 14 bit-depth Roper Coolsnaps 
HQ camera. The field of view is 150um x Ilium and 
contains approximately 20,000 particles. The accuracy 
of feature identification in the x-y plane is approximately 
15nm which corresponds to approximately 0.1 of a pixel. 

Due to the confocal slice having a finite thickness, w, 
particles both above and below the focal plane are imaged 
as being in the plane. This introduces an uncertainty, 
w/2 « 400nm, in the z position of the particle. As a re- 
sult the measured distance between two particles, which 
is the distance projected onto the x — y plane, will be 
less than the true distance. The difference between the 
measured and true distance becomes small when r ^> w 
but can be significant at small r. This introduces a sys- 
tematic error in g(r) which can be corrected [26| 
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FIG. 3. The pair-correlation functions, g(r), for a range of 
packing fractions, <j). The data has been corrected for the 
projection effects discussed in the text. The bottom curve is 
for cj) < 0c, for all other curves (j> > (j> c . The red curve is for 
<f> — (j)* . The position of the first peak remains fixed even 
though the nominal particle diameter, indicated by the dots 
at each packing fraction, varies. 



III. RESULTS AND DISCUSSION 



where g raw is the pair-correlation function measured prior 
to taking into account the effects of projection and r > w. 
We have used simulations to check the validity of this 
approximation. 

The correction shifts peaks outward, increases peak 
heights and decreases the valleys of g(r). It primarily 
affects the first peak of g(r). This shown in Fig. [2j The 
position of the first peak is shifted by 4% and its height 
is enhanced by 10%; the effect on higher order peaks is 
significantly smaller. At large r the primary source of the 
noise in g(r) is from under-sampling configuration space. 
The projection is present in all of our data sets, thus will 
not effect relative measurements as a function of (j). 

The data is processed using a locally devel- 
oped C++ implementation of the Crocker- Grier fea- 
ture identification and tracking algorithm [27] in- 
corporating an existing implementation of the iden- 
tification routine [28j|. The tracking and correla- 
tions software is locally developed and is available 
at |https : //github . com/ tacaswell/tracking[ Our 
tracking and correlation code is an order of magnitude 
faster than equivalent code written in an interpreted lan- 
guage. The implementation of spatial correlation func- 
tions scales as 0(n x r^ ax ) where d is the dimension of 
the correlation and n is the total number of particles. 



We will first discuss some of the general features of 
g(r) and then in the subsections below focus on specific 
aspects of the data. 

Figure [3] shows the pair-correlation function at differ- 
ent the packing fractions, <p. This data has been cor- 
rected for the projection issues as discussed in Section ITTI 
In each case, there is a large first peak at n ~ 0.9/ira. At 
the lowest (j) shown, only the first three peaks are large 
enough to be easily discerned. As <p is increased, peaks 
at larger r grow and become clearly visible. 

The position of the first peak, n, does not shift appre- 
ciably even though the packing fraction varies by a factor 
of 1.7. To emphasize this point, the nominal particle di- 
ameter, <j, is marked on each curve in Figure [3l At all 
4> > C , a > ri, indicating that the particles strongly 
interact with and deform their neighbors. As we vary cr, 
or equivalently (/>, the height of g\ varies but the peak 
postions do not move in location. This can be under- 
stood because the average interparticle spacing in this 
densely packed system is set by the number density, n 
which is held constant by the constraint that the number 
of particles in the rigid sample cell is not allowed to vary. 

To determine the absolute packing fraction for our sys- 
tem, we first determine the number density of our system 
since, at <p c and T = 0, t*i = a. Substituting into equa- 
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FIG. 4. Peak heights versus normalized packing fraction. 
gm — 1 is plotted versus 0/0* for the first four peaks, m — 
1, 4. □ are extracted from data where the volume fraction 
is continuously ramped to higher values. • is extracted from 
data quenched to 0. At all volume fractions the vertical error 
bars are on the order of the marker size. For > 0* the data 
points are tightly spaced. 
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tion (pQ) we find 



Q(j) c 



and 



5 10 
peak # [m] 



5 10 
peak # [m] 



(3) 



We use (j) c = 0.64, which is the density of random close 
packing as measured in simulations and hard sphere ex- 
periments [H, . Our samples range in packing fractions 
from 6 0.55 to 6 = 1.17. 



FIG. 5. a) The pair-correlation function for — 0.980*. The 
raw g(r) is shown in green and the corrected g(r), eq ([2]), is 
shown in black. The inset is a zoom of 8um< r < 12um to 
show the detail of the higher order peaks, b) Peak height 
versus peak number m. Peak height are fit to the decay form 
of eq (|4|), shown in blue, c) Peak location versus m. The peak 
locations are evenly spaced and can be fit to straight line, eq. 
(0, with a slope of 0.85. 



A. Behavior of First peak of g(r): Vestige of the 
jamming transition 

In figure HI we plot the height of the m th peak, g m , of 
the first four peaks in g(r) versus 0. We measured gm in 
two protocols. One is to quench the sample rapidly to a 
desired packing fraction and the other is slowly ramp 0. 
The two protocols yield nearly identical results showing 
that the data is not significantly affected by aging or by 
transient effects. 

The first peak, gi, has non-monotonic behavior. Start- 
ing at small 0, g\ grows with until it reaches a max- 
imum value at = 0* and then decreases as is in- 
creased further. The other peaks, m = 2,3,4, appear to 
reach a plateau above 0*. Within the error bars of our 
measurement, they show no peak. The behavior in q j is 
consistent with experiments [20| and simulations pjl, [30| 



on thermal two-dimensional bi-disperse systems. As in 
those cases, we interpret the peak in g\ as a vestige of 
the T = jamming transition. 

The behavior of g\ is due to the interplay between 
particle overlap caused by thermal motion and particle 
overlap due to geometric frustration, as described in (20j . 
This leads to a peak in g\ as a function of <p as observed. 
The maximum in g\ is necessarily shifted to <j) > (f) c by 
the thermal motion. In our three-dimensional samples, 
0* is shifted above C by a much greater amount than 
was reported in the two-dimensional studies. 
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B. Structure: long-range correlations 

The data in Figure [3] shows no sign of a split sec- 
ond peak. This is in marked contrast to studies of hard 
spheres Q and simulations at T = [H, [3l[. In those 
systems, the second peak splits into two sub-peaks as 
the packing fraction approaches <\> c . This splitting has 
been used as a possible signature to identify colloidal 
glasses [2l|. However, we do not see the second peak 
split at any <\> in our amorphous packings. 

At large (/>, many peaks, corresponding to higher-order 
coordination shells, are visible. The height of these peaks 
decays at large r. Remarkably, as shown in Figure [5j 
near cjf we can identify at least 14 peaks in g(r), the 
smallest of which are fluctuations of less than 1% from 
uniform density as seen in the inset. We are unaware of 
any other experiment that identifies as many peaks in an 
amorphous sample. 

To study the form of g(r) we plot the peak height, g m , 
and location, r m , against peak number m in Figure [5] b) 
and c) respectively. We can fit the decay envelope to the 
Percus-Yevic asymptotic form, 



1 = C- 



exp 



(4) 



where r m is the location of the m th peak and C and £ 
are fitting parameters which depend on <p [32|. At </>* 
C = 2.31 and £ = (—2.6 ± .l)ri, which is a longer cor- 
relation length than seen in hard sphere simulations and 
experiments [33l l34j. 

Beyond the first peak, the inter-peak spacing is very 
uniform: the location of the m th peak, r m , is accurately 
given by 



^1 = 1 + (0.85 ±0.02)(ra-l) 



(5) 



as shown in Figure [5] c. This is also true for the higher- 
order peaks at other values of (j) where they are vis- 
ible. The peak spacing is in quantitative agreement 
with what has been observed in other systems: hard- 
sphere colloids [34|, ball bearings [29|, simulations at 
4>c [S3[ and with experimental measurements of liquid 
noble gases [29|, despite the drastically different inter- 
particle potentials and temperatures. This suggest a fun- 
damental geometric origin for this peak spacing. 

These two results taken together suggest that beyond 
the first peak, g(r) can be approximated by a damped 
sinusoid: 



]{r) 



C 



■ exp 



COS 27T 



ri 



(6) 



This form is roughly consistent with a wide range of an- 
alytic [35|| , simulation [33|, [36| , and experimental [34| re- 
sults and indicative of a dense fluid structure. 
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FIG. 6. a) van Hove correlation functions for r = 1.5s. Curves 
for cj) < (j)* are dashed (- -) and curves for <j) > 0* are solid (- 
). The curves for (j) > 1.090* cannot be distinguished on this 
plot because they lie under the curve at 1.070*. The vertical 
dotted lines indicate the average particle radius, b) The full 
width at 1% max as a function of 6. 
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FIG. 7. Mean squared particle displacement as a function 
of r for different values of (j). The dotted line indicates the 
diffusion in a dilute sample. The symbols and colors match 
those in Figure [6] 



C. Dynamics: Slowing at 0* 

We now ask whether packing fraction, 0*, where g\ has 
a peak, is associated with any change in the dynamics in 
the sample. To quantify the particle mobility we use the 
van Hove correlation function, 



P X (A, t) = 1 ( J2 5 -Xi{t + r)- A) 



(7) 



where Xi(t) is the x component of the zth particle location 
at time t, () is an average over all starting times, and N 
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is a normalization constant such that J dA P X (A, r) = 1. 
Physically P x (A,t) is the probability of a particle mov- 
ing a distance A in the x direction in a time r. We 
can similarly define P y (A,r) for motion in the y direc- 
tion. The mobility is isotropic in the x and y directions, 
which allows us to average these distributions to improve 
statistics at large A. We call the averaged distribution 
P(A,t). We can compute the three-dimensional mean- 
squared displacement as the second moment of P(A,r), 



A 2 (r) = 3 / dA A 2 P(A,r). 



(8) 



For a diffusive system P(A, r) is a Gaussian centered on 
A = and A^jY) oc r a with a = 1. 

There is a qualitative change in particle mobility as the 
packing fraction approaches 0*. This change is evident 
is the shape of P(A, r = 1.5s) shown in Figure [6] Curves 
for < 0* and > 0* are plotted as dashed and solid 
lines respectively. The vertical dashed lines are separated 
by ri and indicate at which point particles have moved 
half the distance to their neighbors, ri/2, in 1.5s. 

At the lowest plotted, 0.930*, the distribution is 
nearly Gaussian and a significant number of particles 
move more than ri/2, indicating substantial rearrange- 
ment of the packing on the 1.5s time scale. We can- 
not take dynamical measurements below 0.930* , because, 
over the course of a few minutes, crystalline regions be- 
gin to nucleate. As is increased, P(A) narrows. Here 
we only report dynamic measurements where the sample 
is stable over the entire course of the experiment. Static 
measurements can be made from a single snapshot and 
can thus be acquired from meta- stable configurations. 
Above 1.070* the curves are nearly indistinguishable and 
lie on top of each other. In this dense state, the particles 
are essentially arrested on the time scales probed. Even 
at these large values of 0, the error in finding particle 
positions is much smaller than the width of the distribu- 
tions measured. 

To quantify the change in particle mobility, we com- 
pute A 2 (t), plotted in Figure [7^. For comparison, the 
expected A 2 (r) for a dilute sample is shown as the dot- 
ted line. At the lowest the particles are nearly diffusive, 
but with the diffusion constant significantly reduced from 
the dilute limit. The narrowing of P(A, r) is reflected in 
the vertical shift of the curves and is accompanied by a 
suppression of the slope. If we assume that A 2 (r) oc r a , 
then the by fitting the slope we can extract a(0), shown 
in Figure 03. There is clearly a dynamical transition in 
the neighborhood of 0* , but it does not have a sharp sig- 
nature. Surprisingly, there is a range of where the par- 
ticles are in contact with their neighbors, i.e., > C , but 
the system is nevertheless nearly diffusive. It is tempting 
to associate the dynamical arrest at 0* with the vestige 
of the jamming transition. However this is problematic 
because the peak in g\ is an unambiguous geometric sig- 
nature and does not display aging effects, whereas the 
dynamical signature is not sharp and depends on the ex- 
perimental timescale used to perform the measurements. 




FIG. 8. The pair-correlation function, g(r), at 0.960* and 
1.060* (on left axis). The difference between these curves, 
Ag, (right axis) is smaller than the experimental resolution 
despite the drastic difference in particle mobility. 



It has been argued that the particle mobility can be 
derived from the pair-correlation function [37j . However, 
if we compare g(r) and a at 0.960* and 1.050*, we see 
that despite an order of magnitude difference in particle 
mobility, the structure is experimentally indistinguish- 
able. This is in agreement with recent work that showed 
that, by using different potentials, configurations with 
very similar g(r) can have very different dynamics (38^ . 
These results are potentially problematic for mode cou- 
pling theory as we have experimentally shown that the 
same static structure can give rise to very different dy- 
namics [39|. 



IV. CONCLUSIONS 

We have demonstrated that there is a peak in g\ as 
a function of in a three-dimensional packing of soft 
pNIPAM colloids undergoing Brownian motion. This is 
a vestige of the T = jamming transition that survives 
at finite temperature and is consistent with previous ex- 
periments in two-dimensional systems [20]. However, in 
contrast with the results in two-dimensions, the packing 
fraction, 0*, where g\ has a peak is significantly higher 
than 0c, the packing fraction where particles would first 
jam at zero temperature. This indicates that, although 
the system is very far from the jamming point, aspects of 
the jamming transition are still observed in the structure 
of the sample. We have also shown that our soft colloid 
fluid will become dynamically arrested when the packing 
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fraction is increased above <p c . This is an example of a 
pressure induced glass transition 14(1-43] . We observe 
that the packing fraction where dynamical arrest occurs 
is close to but not identical with (/>*. 

We have also observed that the pair-correlation func- 
tion for this soft-sphere system has features expected of 
a fluid. For example, the second peak in g(r) is smooth 
with no sign of any splitting that is characteristic of hard- 
sphere systems. Moreover, near (/>*, we observe up to 14 
equally spaced peaks in g(r) whose amplitude decays as 
a function of distance in a fashion that is also consis- 
tent with predictions of liquid structure. As the packing 
fraction is decreased below <fi* , the damping of the peaks 
becomes much more dramatic so that near <p c only three 
peaks are clearly visible. We observe that in Fig. 4, 
gi increases rapidly at low <p until it reaches a peak, but 
then decreases more gradually above </>* . This asymmetry 
around the peak can be related to how the contributions 
to the overlap of particles from thermal motion and from 
pressure vary with packing fraction. The role of temper- 
ature in broadening the first peak rapidly becomes less 
important as the packing fraction is increased. 

Because the pNIPAM colloids are so soft, it is possi- 
ble to have very large particle overlap. This not only 



changes the structure of g(r) but also allows the system 
to be diffusive at packing fractions, </>, that are inacces- 
sible to hard-spheres. To understand thermal systems at 
high densities, such as molecular glass formers, it is im- 
portant to understand the dynamics of soft-sphere pack- 
ings. Further studies are needed on the dynamics in the 
neighborhood of </>*, particularly in comparison to the 
dynamical heterogeneity and correlated motion observed 
in hard-sphere packings near (f) c 8} . This system can also 
be used to measure the density of states to determine if 
the dynamical predictions at T = extend to T > in 
three dimensions. 
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